Effect of cosmic backreaction on the future evolution of an accelerating universe 
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We investigate the effect of backreaction due to inhomogeneities on the evolution of the present 
universe by considering a two-scale model within the Buchert framework. Taking the observed 
present acceleration of the universe as an essential input, we study the effect of inhomogeneities in 
the future evolution. We find that the backreaction from inhomogeneities causes the acceleration 
to slow down in the future for a range of initial configurations and model parameters. The present 
acceleration ensures formation of the cosmic event horizon, and our analysis brings out how the 
effect of the event horizon could further curtail the global acceleration, and even lead in certain 
cases to the emergence of a future decelerating epoch. 
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I. INTRODUCTION 

The present acceleration of the Universe is well estab- 
lished observationally[l], but far from understood theo- 
retically. The simplest possible explanation provided by a 
cosmological constant is endowed with several conceptual 
problems [ J] . Alternative scenarios based on mechanisms 
such as modification of the gravitational theory and in- 
voking extra fields invariably involve making unnatural 
assumptions on the model parameters, though there is 
no dearth of innovative ideas to account for the present 
acceleration [3]. The present state of knowledge offers 
no clear indication on the nature of the big rip that our 
Universe seems to be headed for. 

Observations of the large scale universe confirm the 
existence of matter inhomogeneities up to the scales of 
superclusters of galaxies, a feature that calls for at least 
an, in principle, modification of the cosmological frame- 
work based on the assumption of a globally smooth 
Friedmann-Robertson-Walker (FRW) metric. Taking 
the global average of the Einstein tensor is unlikely to 
lead to the same results as taking the average over all 
the different local metrics and then computing the global 
Einstein tensor for a nonlinear theory such as general 
relativity. This realization has lead to investigation of 
the question of how backreaction originating from den- 
sity inhomogeneities could modify the evolution of the 
universe as described by the background FRW metric at 
large scales. 

In recent times there is an upsurge of interest on study- 
ing the effects of inhomogeneities on the expansion of the 
Universe. The main obstacle to this investigation is the 
difficulty of solving the Einstein equations for an inhomo- 
geneous matter distribution and calculating its effect on 
the evolution of the Universe through tensorial averaging 
techniques. Approaches have been developed to calcu- 
late the effects of inhomogeneous matter distribution on 
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the evolution of the Universe, such as Zalaletdinov's fully 
covariant macroscopic gravity [')]: Buchert's approach of 
averaging the scalar parts of Einstein's equations [(i, 7] 
and the perturbation techniques proposed by Kolb et. al. 
[8] . Based on the framework developed by Buchert it has 
been argued by Rasanen [!)] that backreaction from in- 
homogeneities from the era of structure formation could 
lead to an accelerated expansion of the Universe. The 
Buchert framework from a different perspective devel- 
oped by Wiltshire [10] also leads to an apparent accel- 
eration due to the different lapse of time in underdense 
and overdense regions. 

Backreaction from inhomogeneities provides an inter- 
esting platform for investigating the issue of acceleration 
in the current epoch without invoking additional physics, 
since the effects of backreaction gain strength as the in- 
homogeneities develop into structures around the present 
era. It needs to be mentioned here that the impact of in- 
homogeneities on observables [11, 12] of an overall homo- 
geneous FRW model has been debated in the literature 
[13[. Similar questions have also arisen with regard to 
the magnitude of backreaction modulated by the effect 
of shear between overdense and underdense regions [14[. 
While further investigations with the input of data from 
future observational probes are required for a conclusive 
picture to emerge, recent studies [6-10, 13, 15, 16] have 
provided a strong motivation for exploring further the 
role of inhomogeneities in the evolution of the present 
Universe. 

The backreaction scenario within the Buchert frame- 
work ]6, 15, 16] has been recently studied by us from 
new perspective, viz., the effect of the event horizon on 
cosmological backreaction ]17]. The acceleration of the 
universe leads to a future event horizon from beyond 
which it is not possible for any signal to reach us. The 
currently accelerating epoch dictates the existence of an 
event horizon since the transition from the previously 
matter-dominated decelerating expansion. Since backre- 
action is evaluated from the global distribution of matter 
inhomogeneities, the event horizon demarcates the spa- 
tial regions which are causally connected to us, and hence 
impact the evolution of our part of the Universe. By con- 



sidering the Buchert framework with the exphcit presence 
of an event horizon during the present accelerating era, 
our results [17] indicated the possibility of a transition to 
a future decelerated era. 

In the present paper we perforin a comprehensive anal- 
ysis of various aspects of the above model for the future 
evolution of our universe in the presence of cosmic inho- 
mogeneities. We consider the universe as a global domain 
V which is large enough to have a scale of homogeneity to 
be associated with it. This global domain V is then fur- 
ther partitioned into overdense and underdense regions 
called M and £ respectively, following the approach in 
[16]. In this setup we explore the future evolution of 
the universe and then compare the results to those that 
are obtained when the effect the event horizon is taken 
into account. In this manner we are able to demarcate 
the regions in parameter space for which the effect of 
backreaction from inhomogeneities causes the global ac- 
celeration to first slow down, and then ensures the onset 
of another future decelerating era. 

The paper is organized as follows. In Section II we 
briefly recapitulate the essential details of the Buchert 
framework [6, 7, 15] including the two-scale model of 
inhomogeneities presented in [!(>]. Next, in Section III 
we follow the approach of [Ki] and investigate the fu- 
ture evolution of the universe assuming its present stage 
of global acceleration. We then illustrate the effects of 
the event horizon by first assuming an exponential ex- 
pansion for the universe, and later analyzing numerically 
the evolution equations in Section IV. Subsequently, in 
Section V we perform a quantitative comparison of var- 
ious dynamical features of our model in the presence of 
an explicit event horizon [17] with the future evolution 
as obtained within the standard Buchert framework [16]. 
Finally, we summarize our results and make some con- 
cluding remarks in Section V. 



II. THE BACKREACTION FRAMEWORK 

A. Averaged Einstein equations 

In the framework developed by Buchert [6, 7, 15, 16] for 
the Universe filled with an irrotational fluid of dust, the 
space-time is foliated into flow-orthogonal hypersurfaces 
featuring the line-element 



ds^ = ~dt^ 



9zj 



dX\lX^ 



(1) 



where the proper time t labels the hypersurfaces and 
X* are Gaussian normal coordinates (locating free-falling 
fluid elements or generalized fundamental observers) in 
the hypersurfaces, and g^^ is the full inhomogeneous three 
metric of the hypersurfaces of constant proper time. The 
framework is applicable to perfect fluid matter models. 

For a compact spatial domain V, comoving with the 
fluid, there is one fundamental quantity characterizing it 



and that is its volume. This volume is given by 

\V\g = / d^lg 

Jv 



(2) 



where d/ig = y^(^T^(J^X\J{^~X^dX^dX^dX^. From 
the domain's volume one may define a scale-factor 
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1/3 



(3) 



that encodes the average stretch of all directions of the 
domain. 

Using the Einstein equations, with a pressure-less fluid 
source, we get the following equations [6, 15, 16] 
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Here the average of the scalar quantities on the domain 
V is defined as, 
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and where p, TZ and Hx> denote the local matter density, 
the Ricci-scalar of the three- metric gij , and the domain 
dependent Hubble rate iJ-p = a-D/a-D respectively. The 
kinematical backreaction Q-p is defined as 



Qv 
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(8) 
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where 9 is the local expansion rate and a" — _l/^i^jjl 
is the squared rate of shear. It should be noted that Hxi 
is defined as iJp = 1/3 {6) j,. Q-p encodes the depar- 
ture from homogeneity and for a homogeneous domain 
its value is zero. 

One also has an integrability condition that is neces- 
sary to yield (5) from (4) and that relation reads as 
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From this equation we see that the evolution of the back- 
reaction term, and hence extrinsic curvature inhomo- 
geneities, is coupled to the average intrinsic curvature. 
Unlike the FRW evolution equations where the curva- 
tme term is restricted to an a^^ behaviour, here it is 
more dynamical because it can be any function of a-p ■ 



B. Separation formulae for arbitrary partitions 

The "global" domain V is assumed to be separated into 
subregions J^i , which themselves consist of elementary 



space entities Tg that may be associated with some 
averaging length scale. In mathematical terms 2? = U^J^, 
where Tg = U„ J"^^"^ and J"^^"^ n J"i^^ = for &l\ a ^ /3 
and £ ^ m. The average of the scalar valued function / 
on the domain I) (7) may then be split into the averages 
of / on the subregions Tg in the form, 
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\'D\g, is the volume fraction of the sub- 



where \i = \Ti 
region J-g. The above equation directly provides the ex- 
pression for the separation of the scalar quantities p, TZ 
and Hj) = 1/3 (0)p. However. Q-p, as defined in (8), 
does not split in such a simple manner due to the (0)^ 
term. Instead the correct formula turns out to be 



Q-D = ^\iQi 
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where Qi and H^ are defined in J^£ in the same way as Qd 
and Hxi are defined in V. The shear part (c^) ^ is com- 
pletely absorbed in Qi , whereas the variance of the local 
expansion rates(0^) — (6)2, is partly contained in Qe but 

also generates the extra term 3 J2e^m ^i\n {Hi — H.„^) . 
This is because the part of the variance that is present in 
Qg, namely {d'^) -p — {&) jr only takes into account points 
inside Ti. To restore the variance that comes from com- 
bining points of Fi with others in J>„, the extra term 
containing the averaged Hubble rate emerges. Note here 
that the above formulation of the backreaction holds in 
the case when there is no interaction between the over- 
dense and the underdense subregions. 

Analogous to the scale-factor for the global domain, 
a scale-factor at for each of the subregions J-g can be 
defined such that \'D\g — ^^ \^i\g: a-nd hence. 
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Ea. 
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where Xe 



|J>, 



IP, In is the initial volume fraction 



of the subregion J-"^. If we now twice differentiate this 
equation with respect to the foliation time and use the 
result for dg from (5), we then get the expression that 
relates the acceleration of the global domain to that of 
the sub-domains: 
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Following the simplifying assumption of [lii] (which cap- 
tures the essential physics), we work with only two subre- 
gions. Clubbing those parts of V which consist of initial 
overdensity as Ai (called 'wall'), and those with initial 
underdensity as £ (called 'void'), such that D = AiU £, 
one obtains Hx> = XmHm + ^sHs, with similar expres- 
sions for (/o)p and {TVj-p, and 



— — Xm ^ ^£ 1~ 2Xm^£{Hm ^ HsY 
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III. FUTURE EVOLUTION WITHIN THE 
BUCHERT FRAMEWORK 

We now try to see what happens to the evolution of 
the universe once the present stage of acceleration sets in. 
Note, henceforth, we do not need to necessarily assume 
that the acceleration is due to backreaction [9, 1()]. For 
the purpose of our present analysis, it suffices to consider 
the observed accelerated phase of the universe [4] that 
could occur due to any of a variety of mechanisms [.'>] . 

Since the global domain T) is large enough for a scale 
of homogeneity to be associated with it, one can write, 



Pis = / V^d'X = .f{r)aUt) 
Jv 



(15) 



where /(r) is a function of the FRW comoving radial 
coordinate r. It then follows that 
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dp =: CpO,p^ 



(16) 



and hence the volume average scale-factor av and the 
FRW scale-factor ap are related by av ~ cpap, where 
cp is a constant in time. Thus, Hp « Hv, where Hp is 
the FRW Hubble parameter associated with V. Though 
in general Hv and Hp could differ on even large scales 
[Ki], the above approximation is valid for small metric 
perturbations. 

Following the Buchert framework [6, 16] as discussed 
above, the global domain V is divided into a collection of 
overdense regions A4 ~ UjAi^ , with total volume \Ai\g = 
X); l-^"'!'/ ! and underdense regions £ ~ Uj£^ with total 
volume \£\g = ^ \£''\g- Assuming that the scale-factors 
of the regions £^ and M^ are, respectively, given by a^ = 

cs-t" and umj ~ CMjt^ where a, /3, C£- and cmj are 
constants, one has 



^3 _ ^3 .3a. 
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is a new constant, and similarly 
The volume fraction of the subdomain A4 is 
- — —, which can be rewritten in terms of 



\v\ 



where ( 
for Cm 
given by Aa4 

the corresponding scale factors as Xm ~ °-^ L, , 

therefore find that the global acceleration equation (14) 
becomes 
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Figure 1: Global acceleration parameter plotted vs. time. 
The parameter values used are: (i) a = 0.995, /3 = 0.5, (ii) 
a = 0.984, /3 = 0.5, (iii) a = 1.02, /3 = 0.66, (iv) a = 1.02, 
P = 0.66. (The curves (i) and (iii) correspond to the case 
when an event horizon is included in the analysis in Section 
IV). 



where gMh i^ ^ constant. We obtain numerical solutions 
of the above equation for various parameter values (see 
curves (ii) and (iv) of Fig.l). We later compare these 
with the solutions of the model with an explicit event 
horizon studied in Section IV. 



A. Backreaction and Scalar Curvature 

It is of interest to study separately the behaviour of 
the backreaction term in the Buchert model [(i, !(>]. The 
backreaction Qp is obtained from (4) to be 






(19) 



Note that we are not considering the presence of any 
cosmological constant A as shown in (4). We can as- 
sume that (p)p behaves like the matter energy density, 
i.e. {p)i) = -%-, where Cp is a constant. Now, obser- 
vations tell us that the current matter energy density 
fraction (baryonic and dark matter) is about 27% and 
that of dark energy is about 73%. Assuming the dark 
energy density to be of the order of 10~^^ {GeV) , we 
get pvo ^ 3.699 x lO'^^ {GeVf. Thus, using the values 
for the global acceleration computed numerically, the fu- 
ture evolution of the backreaction term Q-p can also be 
computed (see curves (ii) and (iv) of Fig. 2, where we have 
plotted the backreaction density fraction 51 g = ^'Mp')- 
Once we compute the backreaction it is straightforward 
to calculate the scalar curvature {TVj^y as 



Figure 2: Global backreaction density fraction vs. time. The 
parameter values used are the same as in Fig. 1. 



B. Effective Equation of State 

To point out the analogy with the Friedmann equations 
one may write (4) and (5) as 



a-D 



3^ = -A7:G{pfff + 3pf^j) 
iHl = SnGp, 



eff 



(21) 



where the effective energy density and pressure are de- 
fined as 
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Note that as in (19), here also we are not considering the 
presence of any cosmological constant A. In this sense, 
Qx) and (TZ) p may be combined to some kind of dark fluid 
component that is commonly referred to as X-matter. 
One quantity characterizing this AT-matter is its equation 
of state given by 
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Peff 
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(23) 

which is an effective one due to the fact that backreaction 
and curvature give rise to effective energy density and 
pressure. We again plot this effective equation of state 
by computing its value numerically (see curves (ii) and 
(iv) of Fig. 3), and in the next Section we compare the 
reults with those obtained by considering the effect of an 
event horizon. 



IV. EFFECT OF EVENT HORIZON 



(7e)j, = 16^G(p), 
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l^"J In the previous section we have performed a detailed 

analysis of the future evolution of an accelerating uni- 
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Figure 3: Effective equation of state vs. time. The parameter 
values used are: (i) a = 0.995, /3 = 0.5, (ii) a = 0.984, 
P = 0.5, (iii) a = 1.02, /3 = 0.66, (iv) a = 1.02, f3 = 0.66 



Figure 4: The dimensionless global acceleration parameter 
plotted vs. time (s). The parameter values for curves (i) and 
(ii) are a = 0.995, /3 = 0.5, and for curves (iii) and (iv) are 
a = 1.02, P = 0.66. 



verse by calculating the effect of backreaction from inlio- 
mogeneities on the various dynamical quantities within 
the context of the Buchert framework. As mentioned ear- 
Her, an accelerated expansion dictates the formation of 
an event horizon, and our earlier analysis may be inter- 
preted as a tacit assumption of a horizon scale being set 
by the scale of homogeneity labelled by the global scale- 
factor, i.e. Th ~ a-D- In the present section we introduce 
explicitly such an event horizon. Given that we are un- 
dergoing a stage of acceleration since transition from an 
era of structure formation, our aim here is to explore the 
subsequent evolution of the Universe due to the effects 
of backreaction in presence of the cosmic event horizon. 
Though, in general, spatial and light cone distances and 
corresponding accelerations could be different, as shown 
explicitly in the framework of Lemaitre-Tolman-Bondi 
(LTB) models [18], an approximation for the event hori- 
zon which forms at the onset of acceleration could be 
defined by 
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(24) 



in the same spirit as op = cpap. 

Since an event horizon forms only those regions of T) 
that are within the event horizon are accessible to us. 
Hence, in this case an apparent volume fraction Xmh ^ 
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\MAj, 



is introduced. From (17) it follows that 
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where c^ 



^M\ 



|A^i|g/|7r is a constant. Normalizing 
the total accessible volume in the presence of the event 
horizon, we can write 
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where Af,^ is the apparent volume fraction for the sub- 
domain £. It hence follows that the global acceleration 
equation (14) is now given by 



ay 



CA^.i'^/3(/3-l) 



t2 



1- 



^M 



ff^\ a{a-\) 



f2 



"Mh 



t^f 



1- 




t 



(27) 



In order to obtain the future evolution of the universe 
with backreaction in presence of the event horizon, one 
has to solve the above equation for the scale-factor with 
the event horizon r^ given by (24). 

Though we will eventually obtain numerical solutions 
of the above integro-differential equations, it is instruc- 
tive to first obtain some physical insight of the evolution 
by taking recourse to a simple approximation. To this 
end, let us for the moment model the onset of the present 
acceleration of the Universe by an exponential expansion, 
keeping our analysis close to observations. Specifically, 
we set ax> oa exp (H-pt) in equation (24) only. Using 
Hp = H-p , where Hp is the FRW Hubble parameter 
associated with V, it follows that r^ = Hp^, a constant 
which we substitute in equation (27). With this substi- 
tution, the global acceleration a-p vanishes at times given 
by 

^3/3 = [(3/3- a -1) 



±V(3/S- 



iy + 8a{a-l)] (28) 



The scalefactor of the 'wall' grows as t^ , where 1/2 < (3 < 
2/3. Equation (28) corresponds to real time solutions for 
1 



a> 



(/3 + l) + 2^2/3 (1-/3) . 
Now, let us consider the following two cases separately. 
Case r. a < 1 and /? < 2/3. There exist two real so- 
lutions of (28) corresponding to two values of time when 



the global acceleration vanishes. In Fig. 4, we plot a 
dimensionless global acceleration parameter °° .> with 

° 0,0^5 

time using equation (27). The curve (i) corresponds to 
this case showing that the Universe first enters the epoch 
of acceleration due to backreaction, which subsequently 
slows down and finally vanishes at the onset of another 
decelerating era in the future. 

Case ir. a > 1 and /3 < 2/3. From (28) it follows 
that there is only one real solution (minus sign for the 
square root). This case models the Universe which accel- 
erates due to some other mechanism (not backreaction), 
but subsequently enters an epoch of deceleration due to 
backreaction of inhomogeneities in the presence of the 
event horizon [see curve (iii) of Fig. 4]. 

The plots in Fig. 4 have been done taking the standard 
values of the parameters Vh = H^^ = 4.36xl0^^s, while 
choosing the appropriate range for the parameters a and 
/3, as given in the figure caption. Based on the N-body 
simulation values used in [l(i], we also take Xmho — 0.09. 

We now study the acceleration (27) numerically with- 
out assuming a priori any behaviour for the horizon. 
Keeping with the spirit of our analysis, we assume that 
the Universe has entered the accelerated stage and thus 
a cosmic event horizon has formed. This ensures that 
r^ defined by (24) will be finite valued, enabling us to 
replace the integral equation (24) by 
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Thus, the evolution of the scale-factor is now governed 
by the set of coupled differential equations (27) and (29). 
We numerically integrate these equations by using as an 
'initial condition' the observational constraint go ~ —0.7, 
where go is the current value of the deceleration parame- 
ter (we also used this condition to plot the curves for the 
exponential case), and using the solution for the scale- 
factor plot the global acceleration versus time in Fig. 4 
(curves (ii) and (iv) [thus, all the curves in Fig. 4 are 
set to intersect at the point (io, go)|- The values of the 
other parameters including a and /3 are chosen to be the 
same as in the corresponding curves of the exponential 
case. Curves (ii) and (iv) in Fig. 4 are drawn from the 
numerical solutions where curve (ii) corresponds to the 
case where a < 1 and curve (iv) corresponds to the case 
where a > 1. For a < 1 we see that for the numerical 
case (curve (ii)) the acceleration decreases much more 
sharply than the analytic (exponential) case (curve (i)). 
In contrast, for a > 1, the acceleration which has only 
one transition between acceleration and deceleration in 
the future, looks much more similar for both the (ex- 
ponential and numerical) cases. Therefore, it turns out 
that the exponential expansion for the universe is a much 
better approximation when we have a > 1. The differ- 
ences in the various slopes and also in the scale for the 
dimensionless global acceleration parameter in the two 
cases arise as a result of the approximation of constant 
horizon used in the former, as well as due to the choice 
of the condition go ~ —0.7 used in both the cases. 
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Figure 5: Global scalar curvature density fraction, f?^ plotted 
vs. time. The parameter values for curves (i) and (ii) are 
a = 0.995, /? = 0.5 , and for curves (iii) and (iv) are a — 1.02, 
/3 = 0.66. 




Figure 6: Effective equation of state, Wa^e// plotted vs. time. 
The parameter values for curves (i) and (ii) are a = 0.995, 
/3 = 0.5 , and for curves (iii) and (iv) are a = 1.02, /3 — 0.66. 



We now compute the backreaction Qx, for this model, 
using Eq. (19). Once we compute the backreaction it 
is straightforward to obtain the scalar curvature {TVj-p 
numerically, as was done in the previous section for the 
model without the event horizon. In Fig. 5 we have plot- 
ted the scalar curvature density fraction H,!^ ~ — „2 ■ 
The curves (i) and (iii) correspond to the exponential 
case and curves (ii) and (iv) correspond to the numerical 
case. We observe from Fig. 5 is that the scalar curva- 
ture {TVj-p turns out to be negative when the Universe is 
in an accelerated phase, which is what we expect from 
our knowledge of FRW cosmology. For a < 1 the cur- 
vature increases to a maximum and then decreases, but 
for a > 1 the curvature density keeps on monotonically 
decreasing. Just like we observed in the acceleration plot 
(Fig. 4), here also the assumption of exponential expan- 
sion is a much better approximation when we have a > 1. 



We next compute the effective equation of state using 
Eq.(23). We see that for the duration of positive acceler- 
ation the equation of state remains negative. For a < 1 
the equation of state first reaches a minimum and then 
keeps on rising. This rise is much more rapid for the nu- 
merical case (curve (ii)). For a > 1 the equation of state 
keeps on rising monotonically. In this case it is quite ap- 
parent that the acceleration of the universe is caused by 
something other than backreaction, since the equation of 
state starts from large negative values. Here also we find 
that the assumption of exponential expansion is a much 
more better approximation when a > 1. 



V. DISCUSSIONS 

Let us now compare the nature of acceleration of the 
Universe for the two models described respectively, in 
Sections III and IV. The global acceleration for the two 
models have been plotted in Fig. 1. Here curves (i) and 
(iii) are for the case when an event horizon is included, 
and curves (ii) and (iv) correspond to the case without 
an event horizon. For a < 1 we see that the acceleration 
reaches a much greater value when no event horizon is 
present and that too very quickly, but decreases much 
more gradually than when an event horizon is included. 
This behaviour could be due to the fact that the inclusion 
of the event horizon somehow limits the global volume of 
domain V in such a way that the available volume of the 
underdense region £ is lesser than when an event hori- 
zon is not included. This causes the overdense region 
M to start dominating much earlier and thus causing 
global deceleration much more quickly. When a > 1 the 
acceleration goes asymptotically to zero when an event 
horizon is not included (curve (iv)), which is contrary to 
the case with an event horizon because there the accel- 
eration always becomes negative in the future whether 
a < 1 or a > 1. In fact here also we can see that when 
we include an event horizon in our calculations then there 
is much more rapid deceleration, the reason for which is 
mentioned above. Initially the curves for the two case 
are almost identical, but later on they diverge due to the 
faster deceleration for the case when an event horizon is 
included. 

A similar comparison of the backreaction for the two 
models is presented in Fig. 2. From the expression of 
fig and (8) it can be seen that the backreaction will be 
dominated by the variance of the local expansion rate 9. 
Here we observe that even for the model where an event 
horizon is not included, fl^ is negative for the duration 
over which the acceleration is positive (although fig is 
positive at the extreme edges). For a < 1 the backreac- 
tion density first reaches a minimum and then keeps on 
rising, this rise being much more rapid for the case where 
an event horizon is included (curve (i)). For a > 1 the 
backreaction density keeps on rising monotonically. Here 
curves for both the models are initially very similar but 
later on they diverge. 




0.62 0.64 0.66 0.68 



Figure 7: The range of parameters a and /3 for which future 
deceleration takes place, shown within the respective contours 
for the curves (i) and (iii) corresponding to the case with 
an event horizon, and the curves (ii) and (iv) corresponding 
to the case where an event horizon is not considered. The 
value of Xmo for the curves (i) and (ii) is 0.15 and for curves 
(iii) and (iv) is 0.2. The shaded region corresponds to the 
curve (iv) deinarcating the parameter space for this case when 
acceleration vanishes in finite future time. 



We next consider the effective equation of state w^ ^ff 
(Fig. 3). It remains negative for the entire duration over 
which the acceleration is positive. For a < 1, w^ ^rr first 
reaches a minimum and then keeps on rising. The falling 
to the minimum is much faster for the model where an 
event horizon is not included (curve (ii)) and the increase 
from the minimum being much faster for the case where 
an event horizon is included (curve(i)). For a > 1 the 
curves for the two cases are very similar initially (again 
like backreaction plots) but later on they diverge. 

From our analysis so far it is clear that the accelera- 
tion of the Universe could become becomes negative in 
the future for certain values of the parameters a and /?, 
which represent the growth rates of the scale factors cor- 
responding to the void and wall, respectively. The range 
of values of a and /? for which a future transition to de- 
celeration is possible, is depicted in Fig. 7. We provide 
a contour plot of a versus (3 demarcating the range in 
parameter space (inside of the contours) for which ac- 
celeration vanishes in finite future time. We see that 
initially the curves for the model with an event horizon 
(curves (i) and (iii)) have the same value of a for vari- 
ous values of /3. But later on when the values of a begin 
to change then the curves for the two models (with and 
without event horizon) merge. Since the acceleration has 
no chance of becoming negative when we have a > 1, 
and for the case with no event horizon (the acceleration 
then goes asymptotically to zero), we therefore get the 
maximum limit of a for this case depicted by the line 
a = l. 
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VI. CONCLUSIONS 

To summarize, in this work we have performed a de- 
tailed analysis of the various aspects of the future evolu- 
tion of the presently accelerating universe in the presence 
of matter inhomogeneities. The backreaction of inhomo- 
geneities on the global evolution is calculated within the 
context of the Buchert framework for a two-scale non- 
interacting void- wall model [6, 7, 15, 16]. We first an- 
alyze the future evolution using the Buchert framework 
by computing various dynamical quantities such as the 
global acceleration, strength of backreaction, scalar cur- 
vature and equation of state. Though in this case we do 
not consider explicitly the effect of the event horizon, it 
may be argued that a horizon scale is implicitly set by the 
scale of global homogeneity labelled by the global scale- 
factor. We show that the Buchert framework allows for 
the possibility of the global acceleration vanishing at a 
finite future time, provided that none of the subdomains 
accelerate individually (both a and /3 are less than 1). 

We next consider in detail a model with an explicit 
event horizon, first presented in [17]. The observed 
present acceleration of the universe dictates the occur- 
rence of a future event horizon since the onset of the 
present accelerating era. It may be noted that though 
the event horizon is observer dependent, the symmetry 
of the equation (27) following from (14) esures that our 
analysis would lead to similar conclusions for a 'void' 
centric observer, as it does for a 'wall' centric one. In 
[17] we had shown that the presence of the cosmic event 
horizon causes the acceleration to slow down significantly 
with time. In the present paper in order to understand 
better the underlying physics behind the slowing down 
of the global acceleration, we have explored the nature 
of the global backreaction, scalar curvature and effective 
equation of state. We have the provided a quantitative 
comparison of the evolution of these dynamical quanti- 
ties of this model with the case when an event horizon is 
not included. 

Our analysis confirms that the evolution with backre- 
action from inhomogeneities in the presence of an event 



horizon for some time after the onset of the present ac- 
celerating epoch in the case when one of the individual 
subdomains, i.e., the underdense wall, is accelerating on 
its own, stays close to the standard ACDM model (with 
exponential expansion). On the other hand, in compar- 
ison with the model without an event horizon, we find 
that during the subsequent future evolution the global 
acceleration decreases more quickly when we include an 
event horizon. The reason for this effect is that in the 
latter model an effective reduction of the volume frac- 
tion for the void leads to the overdense region starting to 
dominate much earlier and hence, causes faster deceler- 
ation of the universe. One important difference between 
the acceleration of the two models is that the acceler- 
ation does not vanish in finite time, but in stead goes 
asymptotically to zero for a > \ when no event horizon 
is included. Nonetheless, when a > 1, the curves for ac- 
celeration, backreaction, scalar curvature, and effective 
equation of state for both the cases are very similar and 
only diverge later on. We finally demarcate the region 
in the parameter space of the growth rates of the void 
and the wall, where it is possible to obtain a transition 
to deceleration in the finite future. 

Our results indicate the fascinating possibility of back- 
reaction being responsible for not only the present accel- 
eration as shown in earlier works ]9, Ki], but also leading 
to a transition to another decelerated era in the future. 
Another possibility following from our analysis is of the 
Universe currently accelerating due to a different mech- 
anism [■■)], but with backreaction ]r), 1.5, 10] later causing 
acceleration to slow down. Finally, it may be noted that 
the formalism for obtaining backreaction [G, 7, 1-5, 16] 
used here relies on averaging on constant time hypersur- 
faces. Recently, another procedure based on light-cone 
averaging has been proposed ]19]. Though, the full impli- 
cations on the cosmological scenario of the latter frame- 
work are yet to be deduced, our analysis of including the 
effect of the event horizon introduces certain elements of 
lightcone physics in the study of the role of backreaction 
from inhomogeneities on the global cosmological evolu- 
tion. 



[1] S. Perlmutter, et. al., Nature 391, 51 (1998); A. G. Riess 
et. al., Ap. J. 116, 1009 (1998); M. Hicken et al., ApJ, 
700, 1097 (2009); M. Seikel, D. J. Schwarz, JCAP 02, 
024 (2009). 

[2] S. Weinberg, Rev. Mod. Phys., 61, 1 (1989); T. Padman- 
abhan. Adv. Sci. Lett. 2, 174 (2009). 

[3] V. Sahni, in Papantonopoulos E., ed.. Lecture Notes in 
Physics Vol. 653, The Physics of the Early Universe. 
Springer, Berlin, p. 141 (2004); E. J. Copeland, M. Sami, 
and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 (2006). 

[4] M. Seikel, D. J. Schwarz, J. Cosmol. Astropart. Phys., 
02, 024 (2009) 

[5] R. Zalaletdinov, Gen. Rel. Grav., 24, 1015 (1992); Gen. 
Rel. Grav., 25, 673 (1993). 



[6] T. Buchert, Gen. Rel. Grav., 32, 105 (2000); T. Buchert, 
Gen. Rel. Grav. 33, 1381. (2001). 

[7] T. Buchert, M. Carfora, Phys. Rev. Lett. 90, 031101 
(2003). 

[8] E. W. Kolb, S. Matarrese, A. Notari and A. Riotto, Phys. 
Rev. D 71, 023524 (2005). 

[9] S. Rasanen, J. Cosmol. Astropart. Phys., 0402, 003 
(2004); S. Rasanen, J. Cosmol. Astropart. Phys., 0804, 
026 (2008); S. Rasanen, J. Cosmol. Astropart. Phys., 02, 
on (2009); S. Rasanen, Phys. Rev. D 81, 103512 (2010). 
[10| D. L. Wiltshire, New J. Phys. 9, 377 (2007); D. L. Wilt- 
shire, Phys. Rev. Lett. 99, 251101 (2007). 
[11] A. Ishibashi, R. M. Wald, Classical Quantum Gravity, 
23, 235 (2006). 



[12] A. Paranjape, T. P. Singh, Phys. Rev. D 78, 063522 

(2008); T. P. Singh, arXiv:1105.3450. 
[13] E. W. Kolb, S. Marra, S. Matarrese, Phys. Rev. D, 78, 

103002 (2008). 
[14] M. Mattsson, T. Mattsson, J. Cosmol. Astropart. Phys. 

10, 021 (2010). 
[15] T. Buchert, M. Carfora, Class. Quant. Grav. 25, 195001 

(2008). 
]16] A. Wiegand, T. Buchert, Phys. Rev. D 82, 023523 



(2010). 
[17] N. Bose, A. S. Majumdar, MNRAS: Letters 418: L45-L48 

(2011) 
[18] K. Bolejko, L. Andersson, J. Cosmol. Astropart. Phys., 

10, 003 (2008) 
[19] M. Gasperini et al., J. Cosmol. Astropart. Phys., 07, 008 

(2011) 



